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<^ ■ ABSTRACT 

^ ■ Continuing our series of papers on the three-dimensional (3-D) structures of 

(N ■ and accurate distances to Planetary Nebulae (PNe), we present our study of the 

planetary nebula NGC 6781. For this object we construct a 3-D photoioniza- 
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O ■ tion model and, using the constraints provided by observational data from the 

T^j- ■ literature we determine the detailed 3-D structure of the nebula, the physical 

^ . parameters of the ionizing source and the first precise distance. The procedure 

O ■ consists in simultaneously fitting all the observed emission line morphologies, 

^ . integrated intensities and the 2-D density map from the [SII] line ratios to the 



parameters generated by the model, and in an iterative way obtain the best fit 
for the central star parameters and the distance to NGC 6781, obtaining values 
'^ . of 950 ± 143 pc and 385 L© for the distance and luminosity of the central star 

respectively. Using theoretical evolutionary tracks of intermediate and low mass 



stars, we derive the mass of the central star of NGC 6781 and its progenitor to 
be 0.60 ±O.O3M0 and 1.5 ±O.5M0 respectively. 

Subject headings: Planetary Nebula - Interstellar medium - photoionization mod- 
eling 



1. Introduction 

Planetary nebulae are the end products of the evolution of stars with masses below 
about 8Mq (Pottasch (1984) and Iben & Renzini (1983)). The importance of these objects 
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extends beyond the understanding of how the outer layers of stars end up forming the many 
observed morphologies. Indeed, PNe have been used for many purposes, from understanding 
basic atomic and plasma processes (AUer (1987)) to determining chemical evolution of our 
Galaxy (Maciel & Costa (2003)). Recently PNe are also being used to study galaxies other 
than our own, providing powerful tools to determine distances, kinematics and chemical 
properties of external galaxies (Ciardullo (2003)) and even to trace inter-cluster material 
as in Feldmeier et al. (2004). PNe have the advantage that they are observable out to 
large distances due to their luminous narrow emission lines, especially [OIII]500.7nm, being 
useful also as standard candles through the use of the PN luminosity function (for extensive 
discussion see Jacoby et al. (1999). 

One of the most important problems in observational Galactic PNe research is the diffi- 
culty in determining their distances and three-dimensional structures. Observations always 
produce a 2-D projection of their 3-D structure, and recovering the original structure is not 
trivial. This is also made worse by the fact that only crude distances, usually obtained from 
statistical methods on large samples, can be determined. Large uncertainties are generated 
by the need to assume constancy of one parameter such as the nebular size, (ionized) mass, 
flux etc. so that typical errors in the distances to individual objects are of the order of a 
factor of 3 or more. Very few nebulae have had individual accurate distances determined. 

Historically, PNe have been studied with empirical methods and one dimensional pho- 
toionization models, leading to the above mentioned problems. In contrast to this, our 
technique developed and described in Monteiro, Schwarz, Gruenwald & Heathcote (2004) 
provides precise, self-consistently determined distances, as well as the physical parameters 
for the central star and gaseous nebula, for objects with sufficient observational constraints. 
These objects can provide valuable calibration for existing distance scales as well as self- 
consistently determined physical and chemical quantities. For a detailed description of our 
novel method see our previously published papers in this series and especially the exten- 
sive explanatory appendix in the second paper Monteiro, Schwarz, Gruenwald, Guenthner 
& Heathcote (2005). 

In this work we focus on the PN NGC6781 (RA 19^^ 18™ 28" Dec. +06° 32' 19" 2000.) 
shown in Figure 1. This is a PN whose main structure is observed as a 130" diameter bright 
shell of low ellipticity, double in parts and with fainter lobes emanating at the N-S ends 
where the ring is double and fainter. The ring is brighter in the E-W directions as also 
shown in Figure 1 of Mavromatakis et al. (2001). Their image is available as a fits file on 
the web at http://www.ing.iac.es/ rcorradi and our Figure 1 has been produced using this 
image. 

Mavromatakis et al. (2001) claim a possible faint halo extending out to about 3' by 
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Fig. 1. — Narrowband image of NGC6781 in the light of the [NII]658.4nni hne. Image 
produced from the original of Mavromatakis et al. (2001). 



4' surrounding NGC6781. Note that Corradi et al. (2003) list the object in their paper 
on the search for faint haloes around bright PNe, but do not make any statement about 
having detected a halo according to their criteria, which include the candidate halo having 
to be limb brightened andor detached. The halo of Mavromatakis et al. (2001) is not limb 
brightened and more likely to be the result of scattered light in the instrument. 

For a more detailed discussion of the morphology of NGC 6781 based on a set of narrow 
band images, as well as density, temperature, and extinction maps see Mavromatakis et al. 
(2001). 

An estimate for the luminosity, temperature, and mass of the CS of NGC 6781 based 
on statistical analysis is 127L0, 105 kK, and 0.6 Mq, respectively, taken from Stanghellini, 
Villaver, Manchado & Guerrero (2002). The average distance of the 12 literature values 
we found in Acker et al. (1992) is 946 pc with the individual values ranging from 500 to 
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1600 pc, a factor of more than three. 
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Fig. 2. — The graph shows the Spectral Energy Distribution (SED) of NGC6781 over more 
than five decades of wavelength interval. All values were obtained from the literature. The 
dotted curve is a combined blackbody and continuum radiation fit to the data. 

The Spectral Energy Distribution (SED) AF(A) of NGC6781 shows a broad blackbody- 
like spectrum between about 1 and 100 /im, a radio tail out to 12 cm, and a rise toward the 
blue which probably comes from the hot central star. Integrating the F(A) curve we obtain 
a luminosity of L = 166d^(kpc) Lq. The luminosity of NGC6781 is therefore 150 L© at 
its distance of 950 pc. Applying the usual correction according to Myers et al. (1987) we 
obtain L = 225L0. 

Here we present our own modeling results using observational data published in the 
literature as constraints and derive the 3-D structure, chemical abundances, CS properties, 
and distance for NGC6781 in a self-consistent manner. 

In 52 we discuss the observational data used as constraints for the models. In 53 we 
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present the model results generated by the 3-D photoionization code, and we discuss the 
derived quantities. In §4 we give our overall conclusions and discuss possible discrepancies 
with other determinations of parameters for this nebula. 



2. Observations 

The observational data that we used to constrain our model for NGC6781 were taken 
from the literature. Mavromatakis et al. (2001) presents narrow band imaging of the nebula 
in the most important lines such as Ha ,H/9, [OIII] and others. He also provides [SII] narrow 
band images to obtain a spatially resolved density map, as well as the HaH/? extinction map 
of the nebula. We used their 2-D density map to infer the 3-D density structure to be used 
as an input for the photoionization model. 

Liu, Liu, Luo & Barlow (2004) present deep optical spectra of medium resolution for 
NGC 6781 and 11 other planetary nebulae. The observations were carried out with a long-slit 
spectrograph covering from 360 to 800 nm and include all important emission lines. Their 
observations are particularly interesting because the objects were scanned with the long-slit 
across the nebular surface by driving the telescope differentially in right ascension. These 
observations then yield average spectra for the whole nebula, which are in principle more 
precise than single slit observations, and allow us to produce emission line maps of the object. 
For details of the observations and their reduction procedure as well as the full tables of line 
fluxes, see Liu, Liu, Luo & Barlow (2004). The most important lines used as constraints for 
the model are listed in our Table 1 together with the corresponding values from our model. 

We also used the Ha-I- [NH] image from Mavromatakis et al. (2001) to determine the 
size of NGC 6781 thus using it as one of the constraints for the distance obtained in our 
model calculations. 



3. Photoionization Models for NGC 6781 

The photoionization code we used for the study of NGC 6781 is the Mocassin code 
described in full detail in Ercolano et al. (2003). This code allows for the same possibilities 
as the code used previously in Monteiro, Morisset, Gruenwald & Viegas (2000) but is more 
sophisticated in that the diffuse radiation is fully taken into account in an efficient manner. 
The previous code also had this ability but the associated increase in computational time 
was prohibitive. 

The basic procedure adopted to study NGC 6781 is the same, independent of the code 
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and has been fully described in the two previous papers in this series (Monteiro, Schwarz, 
Gruenwald & Heathcote (2004) and (especially the appendix of) Monteiro, Schwarz, Gruen- 
wald, Guenthner & Heathcote (2005)). In short, we gather as much observational material 
as possible and use it to constrain our model with many data simultaneously. Of partic- 
ular interest are the total line fluxes and line images corrected for reddening, and the line 
diagnostic ratio maps. The structure adopted for the nebulae is defined based on the ob- 
served density map, when available, or density profiles from single slit observations plus the 
observed projected morphology in several emission lines. 

In the case of NGC6781 in particular, the initial structure was based on density maps 
published by Mavromatakis et al. (2001). It is clear from the images and density map (from 
the [SII] doublet ratio) presented here, that the density is lower in the central region than in 
the main bright ring. This indicates that the structure must have lower density material in 
the line of sight of those regions and therefore the best structure to reproduce the observed 
projected morphology is an open ended structure or hour-glass shape, which we therefore 
adopt for this object. It is also clear from the images in many narrow band filters that the 
material is highly clumpy, so to reproduce this we include random density fluctuations in 
the adopted structure. The final adopted structure in its best fitting orientation on the sky 
is presented in Figure 3. 



4. Model Results 

We present here the main results obtained from the photoionization model constrained 
by the observational data. The integrated fluxes for 12 emission lines are given in Table 1, 
as well as the fitted abundances and ionizing star parameters. 

The model fitting procedure which uses the model image size fitted to the observed one 
for the line [NII]658.4nm, as well as the absolute Hj3 flux, and the integrated fluxes of all 
other lines gives a distance of 950±145 pc for NGC 6781. The error on this distance has been 
computed in the same way as in our previous papers in this series. 

In Fig. 4 we show the projected images obtained from the fitted model. Notice that all 
major morphological features of the object are well reproduced as well as the general ioniza- 
tion stratification in the different emission lines. In particular the images for [OIII] 500.7 and 
Hell 468.6 show very good agreement with those obtained by Mavromatakis et al. (2001) 
as well as the more common [Nil] 658.4 nm and Ha. 

Figure 5 shows the final model image for [Nil] 658.4 nm plotted in contours over the 
observed image by Mavromatakis et al. (2001). Again the good agreement of the apparent 



-7- 



size, obtained from the size of the fitted model grid and our determined distance, is evident. 

As with objects studied in previous works, the mass of the ionizing star, as well as its 
progenitor and age are determined from theoretical cooling tracks. Here we have used the 
cooling tracks of Vassiliadis & Wood (1994) because their grids present a good (= close) 
sampling of progenitor masses in the 1 to SMq range. In Figure 6 we show the position of 
NGC6781, as well as the other objects that we have studied previously with this method, 
along with the theoretical cooling tracks. From this we obtain the mass of the central star of 
and its progenitor to be 0.60 ± 0.03 M© and 1.5 ± 0.5 M© respectively. Figure 6 also shows 
the position of these same objects as determined by different authors by distinct techniques. 
All PNe central stars are well evolved on their cooling track except NGC6369. 



5. Discussion and conclusions 

As with all our other PNe central star temperatures determined by this 3-D structure 
method we generally find higher temperatures than those computed from oxygen lines and 
nearer those from Hell lines. Table 2 shows a comparison of parameters from the literature 
and those determined by us in this and previous papers. 

Note that Phillips (2003) found that bipolar PNe have average T(Hell) = 138k; ellipti- 
cals 92 k; round 81k whole sample 87 k. This was also found earlier by Corradi & Schwarz 
(1995) who determined T(bipolars) = 142 k; irregulars 99 k; ellipticals 76 k; and unresolved 
63 k. This general trend is always observed, and our high temperature for bipolar NGC 6781 
is no exception. 

The central star properties of all PNe that we have applied our method to are shown in 
the HR diagram of Figure 6. Also shown are the values determined by other methods, taken 
from the literature. There are large differences between our values and those previously pub- 
lished for NGC 6781 and NGC 6369 and in both cases the central star luminosity determined 
by our method was higher than the literature value. This is because it was assumed that 
these nebulae are radiation bound but we have shown them to be matter bound as they 
lose up to 70% of their UV radiation to space, resulting in an underestimation of both the 
luminosity and temperature. The blue rise in the SED for NGC 6781 also confirms that blue 
radiation is escaping from the nebula. Note that the central star luminosity from our model, 
L = 385 / Lq is larger than the luminosity derived from the observed SED L = 225 /Lq by a 
factor of 1.7, confirming that the nebula is matter bound and that a significant fraction of 
the stellar UV flux escapes from the object. 

We claim that ours are the most accurate luminosities and temperatures that have been 



determined for these stars to date. Interestingly our values tend to bring the core masses 
closer to 0.6 Mq which is also the peak value of the narrow mass distribution of white dwarfs. 

We acknowledge the continuing support of NOAO's Science Fund to HM and the gen- 
erous hospitality of the Nordic Optical Telescope on La Palma. 
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Fig. 3. — The image shows an isodensity surface of the structure adopted for NGC6781 in 
its final fitted spatial orientation indicated in the upper left hand corner. 
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Table 1: Observed and model line fluxes and model central star parameters for NGC6781. 





Observed 


Model 


Rel. Error 


T. (K) 




123kK 




LjLr, 


- 


385 


- 


log {g*) 


- 


7.0 


- 


Distance (pc) 


500-1600 


950 


±0.143 


Density 


100-1400 


100-1400 


- 


He/H 


- 


1.25 X 10-1 


- 


C/H 


- 


5.95 X 10-^ 


- 


N/H 


- 


9.75 X 10-^ 


- 


0/H 


- 


3.50 X 10-^ 


- 


Ne/H 


- 


7.10 X 10-^ 


- 


Ar/H 


- 


2.26 X 10-6 


- 


S/H 


- 


0.28 X 10-5 


- 


\og{H(3) 


-9.80 


-9.84 


-0.09 


Nelll] 386. S'^ 


1.09 


1.06 


-0.03 


OIII] 436.3 


0.05 


0.07 


0.28 


Hell 468.6 


0.09 


0.09 


0.02 


OIII] 500.7 


8.23 


7.10 


-0.14 


Nil] 575.5 


0.07 


0.07 


0.07 


Hel 587.6 


0.16 


0.17 


0.09 


01] 630.2 


0.33 


0.38 


0.15 


Ha656.3 


2.65 


2.86 


0.08 


Nil] 658.4 


3.96 


3.88 


-0.02 


;SII] 671.7 


0.25 


0.25 


-0.01 


SII] 673.1 


0.21 


0.21 


0.00 



a) Value obtained by Kaler, Shaw & Browning (1997) 
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Fig. 4. — Images obtained from the projection at the observed angle of the data cubes of 
emissivities computed by the photoionization code for the most important emission hues. 
Compare these model images with the observed images in Mavromatakis et al. (2001). 
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Fig. 5. — Image comparing the observed [NII]658.4nm narrow band image with the contours 
of the equivalent image of the fitted modeL Note the similarity between the observed and 
modeled images. 
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Fig. 6.— HR diagram for NGC 6781, NGC 3132 (Monteiro, Morisset, Gruenwald & Viegas 
(2000)), NGC 6369 (Monteiro, Schwarz, Gruenwald & Heathcote (2004)), MZ 1 (Monteiro, 
Schwarz, Gruenwald, Guenthner & Heathcote (2005)), all PNe that had their central star 
properties determined by our method. Also plotted are the literature values for comparison, 
a) Stanghellini, Corradi & Schwarz (1993); b) Stanghellini, Villaver, Manchado & Guerrero 
(2002); c) . The evolutionary tracks are from Vassiliadis & Wood (1994); they are similar 
to the Blocker (1995) models but take metalicities etc. into account. . 
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Table 2: Comparison of parameters found in the literature and those determined by our 3-D 
method. 



Object 


Our T 


Hell/HI 


[OIII]/[OII] 


Our L(Lo) 


Our d(pc) 


Hb5 


230^^ 


131^ 




6000^* 


1400" 


Mzl 


120^= 


139^^ 


- 


164^= 


1050^ 


NGC 3132 


90^= 


80^= 


36^ 


150^ 


930^ 


NGC 6369 


91/ 


122^ 


605 


8100^ 


1550^ 


NGC 6781 


123'^ 


126^ 


685 


385'' 


950'^ 



a) Preliminary results from Rice, Schwarz & Monteiro (2004) b) Kaler &: Jacoby (1989) c) Phillips 
(2003) d) Monteiro, Schwarz, Gruenwald, Guenthner & Heathcote (2005) e) Monteiro, Morisset, 
Gruenwald &: Viegas (2000) f) Monteiro, Schwarz, Gruenwald k. Heathcote (2004) g) Gurzadyan 
(1997) h) This paper. 



